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Abstract 

We analyse new contributions to the theoretical input in heavy quark sum rules and we 
show that the general theory of singularities of perturbation theory amplitudes yields 
the method to handle these specific features. In particular we study the inclusion of 
heavy quark radiation by light quarks at 0(a%) and of non-symmetric correlators at 
0(a.s). Closely related, we also propose a solution to the construction of moments of 
the spectral densities at 0{a%) where the presence of massless contributions invalidates 
the standard approach. We circumvent this problem through a new definition of the 
moments, providing an infrared safe and consistent procedure. 
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1 Introduction 



Sum rules analyses have extensively exploited the relation between the correlator of the 
quark electromagnetic currents and the cross section of e + e~ — > hadrons under the assump- 
tion of quark-hadron duality, to extract fundamental information of hadron systems. The 
two-point function containing the QCD dynamics of the produced quarks is built from the 
sum of the electromagnetic vector currents associated to each flavour: 

KaM = ^/rf 4 xe^E^e g ,(0|T(g(x) 7 ^(x)) (7(0) 7 V(0) ) | ) 

<?,<?' 

= (-rp 2 + //)n w (p 2 ), (i) 

where q and q' stand for heavy or light quarks, indistinctly, with electric charges e q and e q >. 
Here we find two types of correlators: the symmetric ones, both electromagnetic currents 
corresponding to the same flavour, and non-symmetric correlators, where q ^ q' . Strictly, 
the latter are needed to fully describe the electromagnetic production of hadrons, even in the 
case where a definite flavour type of hadrons is isolated in the final state. Sum rules analyses 
applied to heavy quark production are written down in terms of the symmetric correlator 
built from the vector current jn(x) = e<g Q(x) 7 M Q{x) of the heavy quark Q, and the effects 
of the non-symmetric correlators are never considered. The reason is that they begin to 
contribute beyond 0(a 2 s ) in QCD perturbation theory (see Fig. 0(a)), which means one order 
beyond the actual knowledge of the (symmetric) heavy quark correlator Hqq- The study of 
such new effects in QQ production will be mandatory if 0(a^) accuracy is reached in the 
future. However, already at 0(a 2 s ) the production of heavy quarks QQ receives contributions 
which have neither been accounted for in the theoretical input of heavy quark sum rules. 
These arise from heavy quark discontinuities of symmetric correlators built from quarks such 
that m q < rriQ, as the cut shown in Fig. 0(b), representing the production of heavy hadrons 
radiated off a pair of lighter quarks. 

Finally, Groote and Pivovarov have recently pointed out [[[], |2[] that, at 0(a^), a three- 
gluon intermediate state (see Fig. 0) contributes to the n Q g correlator. As these authors 
have shown, this massless intermediate state invalidates the usual definition of the moments 
M n , 

1 / A \ n 



(2) 



for n > 4, when they become singular. Consequently the use of heavy quark sum rules at 
0{a1) is debatable. 

All the features we have just quoted arise as a consequence of the interplay between the 
implementation of quark-hadron duality and the proper definition of the observables in the 
case of heavy quarks QCD sum rules. The correlation between the perturbative input and 
the observable information on the experimental side requires a careful matching that cannot 
be fully achieved. Accordingly the introduced incertitudes should be estimated and included 
in the errors of the parameters determined through this method. 
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Figure 1: Examples of perturbative non-heavy quark current correlators at 0(a^) (a) and 
0(a 2 ) (b) that contribute to the production of QQ states. 

Here we discuss the aspects pointed out above and their consequences in the methodology 
of extracting information from QCD sum rules. The aim of this work is to provide a consistent 
procedure to implement the perturbative input in the theoretical side of the heavy quark 
sum rules. Our proposal relies in a careful application of the general theory of singularities of 
perturbation theory. The crucial point will be to isolate all the cuts related to QQ production 
from the general vector two-point function (HD in order to construct a modified correlator 
containing only contributions to heavy quark production. 

In Section 2 we recall the theory of singularities of perturbative amplitudes. The relation 
between the phenomenological and the theoretical input in the QCD sum rules is discussed 
in Section 3. Hence Sections 4 and 5 collect the implementation of our proposal to include 
heavy quark radiation off light quarks and to exclude massless singularities, respectively. We 
will comment on the uncertainties related with our method too. In Section 6 we emphasize 
our conclusions. 

2 Analyticity of U had (s) 

As it is well known two-point functions are analytic except for singularities at simple poles 
or branch cuts, the latter being originated by normal thresholds of production of internal on- 
shell states. Assuming that the absorptive part of Uh a d{p 2 ) starts at some point sq, vanishing 
below this point, the correlator satisfies the dispersion relation [0] [] : 

^ / 9n ^ / 9n ^ / n p 2 f°° ds Imllhad(s) , \ 

n had ( P 2 ) = n had (p 2 ) - u had o = y - \ p-^ . 3 

Ti J s s s — p A — ie 

The absorptive part Imlldd^s) is a physical observable, as it is proportional to the total 
hadron production cross section by a vector current J 1 ' = J2 q jq- Being QCD the underlying 

1 Sometimes the Adler function defined as dH(p 2 )/d\np 2 , to get rid of the subtraction constant, is used. 
The choice of the regularization prescription is not relevant for our discussion here. 
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theory of strong interactions, the quark-hadron duality hypothesis allows us to identify, 
inclusively, the states in terms of observable hadrons with the partonic intermediate states. 
Hence the optical theorem tells us that the total absorptive part is the sum of the absorptive 
parts corresponding to different intermediate partonic states: 



i / E^(0|^|n)(n|jt| ) = J2 

05 J n n 



lmU had (s) = -— I £cLM0| J"|n)(n| 4\°) = E ImII n( s ) > ( 4 ) 

where the phase space integration has been explicitly stated 0. A similar separation between 
contributions of different final hadron states in the perturbative evaluation of the two-point 
correlator, Eq. (|l|), would allow us to keep only the desired heavy quark cuts in the symmetric 
and non-symmetric correlators. Although Cutkosky rules provide a method to isolate cuts 
corresponding to different intermediate states at the perturbative level, some care is needed 
in their application. 

The study of analytic properties of perturbation theory amplitudes shows that their sin- 
gularities are isolated and, therefore, we can discuss each singularity of a perturbative am- 
plitude by itself ||. As a consequence, any one-variable dependent amplitude F(z) satisfies 
a dispersion relation from Cauchy's theorem given by : 

F{2) = _L =E r^l^V)L (5) 

v ' 2m J zf-z V J*n 2vri z' - z ' V ' 

where [-F(z)] n is the discontinuity across a branch cut which starts at the point z n and it 
is associated to a definite intermediate state. For the general two-point function in Eq. (|J), 
which depends on the total momentum squared p 2 , we would have 

rWp 2 ) = E — r — [n(g)]w , (6) 

~ 2ni Js n s s — p 2 — ie ' 

where now [II(s)] n provides the sum of all the cut diagrams associated to a definite inter- 
mediate state labeled n, (n = qq,q'q',ggg,qqq'q', . . .). In the perturbative calculation, every 
discontinuity contributing to [n(s)] can be associated to a "reduced" Feynman diagram 
obtained by contracting internal off-shell propagators to a point and leaving internal on- 
shell lines untouched. Its contribution is written down following the Cutkosky rules for the 
graph. However the discontinuity across a specified cut in a single diagram needs not to 
be a pure real function in the physical region. Hence the separation between the imaginary 
parts coming from different final states, as stated in Eq. (|4|), does not seem to apply for 
individual diagrams. But from Eqs. (^) and we can conclude that [n(s)] n = 2ilmn n (s), 

2 We use dR n — (27r) 4 S 4 (q— X)"=i Pi) Y[i=i dpi, where q is the current four-momentum and dpi = (2nW2E- • 
The — l/(6s) factor in Eq. (|J) originates from Uhad — ~9^u n^ d /(3s) and the (1/2) factor from the unitarity 
relation. 

3 This expression also gives the residue Ri of a pole at z = z% if we interpret the discontinuity as \F(z) ] n = 
2TriRiS(z — Zj). 
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meaning that only the sum of all cuts corresponding to a defined intermediate state provides 
the physical observable, i.e. Imll n (s). Evidently, this holds at any perturbative order in a s , 
and gives a prescription to isolate contributions to different quark intermediate states in the 
hadron two-point function. This assertion might seem obvious but it is not : A QQ cut on 
the right-hand fermion loop in Fig. |2|(a) does not provide, by itself, a pure real contribution. 
Only when both QQ cuts, on the left-hand and right-hand fermion loops of Fig. |2](a), are 
added we get a term contributing to the physical observable Im U n= QQ . 

This last example also shows that some subsets of discontinuities of the same intermedi- 
ate state give already real functions prior to the summation of all contributions at a fixed 
perturbative order. This is the case for the set of cuts coming from a symmetric correlator, 
and for the set arising from a non-symmetric correlator with currents jjf, together with 
its conjugate. This is easily seen if we rewrite the absorptive part corresponding to the state 
n, ImII n , as a sum of terms arising from symmetric and from non-symmetric correlators: 



Iml[ n (s) = - — / (IB,, 



E<°l#l»>W4.l > 



+ E ((0|^|n)(n|4,j0) + (0|jf |n>(n|jtj > 



•(7) 



The first term in the r.h.s. of Eq. (|7|) represents the absorptive contribution from symmetric 
correlators, and the perturbative expansion of each one, following Cutkosky rules, is clearly 
real. In the case of interest, n = [QQ] 0, this term contains the usual heavy quark spectral 
density built from heavy quark currents, IIqq-, and [QQ] production through light quark 
currents correlators. The second and third terms in Eq. (|7|) are conjugate to each other, so 
their sum also gives a pure real number. In terms of diagrams, this means that to extract 
the desired absorptive part from non-symmetric correlators we need to add to the cut of a 
diagram the corresponding one in the conjugated diagram (see Fig. |l](a); the discontinuity 
obtained from the same diagram with quark q and quark Q lines interchanged should be 
added up to get a real contribution). 



3 Phenomenology vs theoretical input in heavy quark 
sum rules 

The analysis above shows that a clear control can be enforced on the perturbative side 
of the sum rules in order to include or exclude specific contributions. However while there 
is no doubt about the observable that provides Imllhad oc c(e + e~ — > hadrons) when an 
exclusive hadron sector (like, for example, heavy quark production) is specified, it is clear that 

4 Brackcts [QQ] are short for any hadron state containing at least a QQ pair and, possibly, light quarks 
and gluons too. 
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the matching between the perturbative and the phenomenological side includes incertitudes 
related with the content and definition of the final state. 

Heavy quark sum rules || have been successful in providing information on the heavy 
quark parameters. In short they make use of global quark-hadron duality that translates 
into the ansatz on the vector correlator ILqq^s) : 

/ ds mil ~ / ds mil + ( 8 ) 

Jso S S 

where ImTtfj^ls) on the l.h.s. gives the phenomenological information on heavy quark 
production and it is related with the cross-section of vector current production of hadrons 
containing Q-flavoured states. On the r.h.s. Im IFj^ (s) is the QCD perturbative contribu- 
tion to the correlator, and in the lower limit of integration M is usually taken as the pole 
mass of the heavy quark. Finally the dots on the r.h.s. are short for non-perturbative (the 
gluon condensate essentially) contributions and possible Coulomb-like bound states coming 
from non-relativistic resummations in n?*f~ below threshold. These last two features are not 
relevant for the discussion of this paper and have to be implemented on our results without 
modification. 

To a definite perturbative order in a s , Imnj^(s) includes all the absorptive contribu- 
tions to the correlator that provide [QCj] production. Notice that this is not the same that 
the absorptive QQ contribution of the heavy-quark current correlator n^g, as it is usually 
assumed. The total experimental cross section cr(e + e~ — > hadrons) can be split into two 
disjoint quantities : the cross section for producing hadrons with Q-flavoured states, and 
the production of hadrons with no Q-flavoured components. If the experimental set up was 
accurate enough to classify events into one of these two clusters, the first class would be the 
required ingredient for the phenomenological part of the heavy quark sum rule. However this 
separation, implemented in the theoretical side within perturbative QCD, is rather involved. 
Up to 0(a 2 s ) there has not been any doubt, in the literature, that contributions to this side 
arise wholly from QQ cuts in the heavy quark correlator Hqq- The physical picture behind 
this assertion relies in the assumption of factorization between hard and soft regions in the 
quark production process and subsequent hadronization. The hard region described with 
perturbative QCD entails the production of the pair of heavy quarks, and the soft part of the 
interaction is responsible for the observed final hadron content. Although possible, annihila- 
tion of the partonic state QQ due to the later interaction is very unlikely, as jets arising from 
the short distance interaction fly apart before long-distance effects become essential. Conse- 
quently, each jet hadronizes to a content of Q-flavoured states with unit probability. As local 
duality is implicitly invoked, this picture is assumed to hold at sufficient high energies; hence 
perturbative corrections to the hard part are successively included through the heavy quark 
currents correlator. We claim, though, that similar QQ cuts are present in non-symmetric 
correlators, starting at 0(a^), as the one shown in Fig. |l](a), where the left hand part of the 
cut diagram is a genuine production of QQ states triggered by virtual light quarks. If the 
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use of heavy quark sum rules up to this order is considered, the inclusion of these terms of 
the correlator of a heavy and a light quark currents should be taken into account. According 
to our conclusion in the last Section, once the discontinuity provided by Fig. |l|(a) is known, 
it has to be added to Im njj^ (s) . 

Other extra QQ cuts, i.e. not contained in IIqq, arise even at 0(a 2 ) as the diagram of 
Fig. 0(b). In this case the QQ pair is produced through the splitting of a hard gluon radiated 
off a pair of light quarks. Whether this cut should be accounted for or not in the theoretical 
side depends crucially on which is the content and the configuration of the reconstructed final 
state in the experimental data, as the physical picture outlined above for pure QQ cuts does 
not apply so clearly for QQqq discontinuities. We will come back to this point at the end of 
Section 4. In addition a discussion about other possible contributing cuts should arise. The 
case of the three-gluon discontinuity is postponed to Section 5. 

In the following we will discuss, in turn, the inclusion of heavy quark radiation by light 
quarks and the infrared massless discontinuities noticed by Groote and Pivovarov. We will 
provide specific solutions along the lines put forward in Sections 2 and 3. 



4 Heavy quark radiation 

Starting at 0(a 2 ), symmetric correlators built from light quark currents include four 
fermion cuts with a heavy quark pair radiated off the light quarks as shown in Fig. 0(b) (two 
additional diagrams, one with the two gluons attached to the lower light fermion line, and 
the other with one gluon attached to each light fermion line, should be considered too). The 
sum of all these four fermion absorptive parts in the three-loop diagrams with massless light 
quarks currents has been calculated in Ref . || , and can be cast into the following form f] : 

a A 2 r s ds' 



12-ImlW*) = R qim = N C (J2 Qi) Cs M) J - R(s') F(s'/s) , (9) 

i=u,d,s 



with 



F(x) = - {(1 + x) 2 ln 2 a; + (3 + 4x + 3a; 2 ) lna; + 5(1 - x 2 ) 
6 I 

r 7T 2 nl 

-4(1+ a;) 2 [Li 2 (-x) + ln(l + x)lnx + — J J . (10) 

The function F(s'/s) gives the rate for the decay of a vector boson of mass y/s into a 
vector boson of mass y/s' plus a pair of massless fermions (qq). The spectral density R(s) = 
(3(3 — (3 2 )/2 (at lowest order) is the normalized cross section for the production of a pair 



of fermions with unit charge through a vector boson; here (3 = yl — AM 2 / s is the velocity 
of the produced heavy quarks. The integral can be solved analytically in this case and the 

5 Notice that our definition of R q qQQ differs from the one in Ref. || . 
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result is found in Ref. 0. Note that the heavy quark pair is created in a colour octet state, 
and the factor 



N c \2 2 J \2 2 J 3 

retains this colour structure. It is interesting to compare the contribution from R q gQQ with the 
0(a 2 s ) contributions to Rqq (i.e. to the spectral density of the heavy quark correlator). Note 
that in the high energy limit there is no difference between the diagram shown in Fig. |I](b) 
and the same one with Q and q lines interchanged or with q = Q, both of them being included 
in U.QQ- Differences arise because the heavy quark currents correlator, Hqq, also accounts 
for two heavy quark cuts where the internal (light or heavy) quark loop represents a virtual 
correction to the electromagnetic current. 

We have written Eq. (§) in terms of a general R(s) function in the integrand because 
it allow us to introduce in a straightforward way final state interactions between the heavy 
quark pair. In particular, we know that close to threshold the Coulomb interaction between 
the heavy quark pair dominates the dynamics. Resummation of leading terms ~ (a s / ' j3) n 
becomes mandatory, and gives rise to the well known Sommerfeld factor multiplying the cross 
section: 

R thr (s) = R(s) x C ™f . (11) 

1 — exp(— C7ia s /p) 

The colour factor C appears in the Coulomb QCD potential and its value depends on the 
relative colour state of the quark pair. For singlet states C = Cp, and the potential is attrac- 
tive, increasing the cross section at threshold. This is the case of heavy quark production in 
e + e _ collisions. However, in our case the heavy quark pair is produced through the splitting 
of a gluon. The Coulomb potential becomes repulsive between quarks in a colour octet state, 
C = Cp — Ca/2 = —1/(2N C ), and the Sommerfeld factor at low velocities then reads 



-7ra s /6/3 ira s 



1 — exp(ira s /Q(3) 6(3 

causing the cross section to decrease near threshold even faster than j3, the phase-space 
velocity in R(s). The production of heavy quarks radiated off massless quarks through a 
virtual gluon is then very much suppressed in the threshold region. However, as mentioned 
above, high energy quark lines can be considered massless and the contribution from this 
diagram is numerically equal to the same one with Q and q lines interchanged. 

The inclusion in ImII??^.(s) of four-fermion cuts coming from light-quark correlators is 
possible because we have shown in Section 2 how to discern and extract these pieces. As 
discussed before, the procedure depends crucially on the definition of the observable infor- 
mation input in the sum rule, and consistence between the theoretical and phenomenological 
parts is required. Let us come back to the discussion of Section 3. There it was argued 
why perturbative QQ cuts are thought to reproduce the phenomenology of two jet events. 
Notice that, in heavy quark radiation from light quarks, the signature of the event is likely 
to be a 3-jet configuration where one of the jets is generated from a gluon. If heavy flavour 
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components are to be found in this jet, the diagram of Fig. |l|(b) would certainly be needed 
to account for these events in the theoretical side. However the heavy partons in this jet are 
not as energetic as in a pure QQ production and, consequently, the proposed factorization 
between long and short distance effects may not longer apply, allowing for an interference 
between both regimes. In this case we cannot argue that these kind of cut diagrams would 
result in a final state with Q-flavoured hadrons with unit probability, although we may im- 
pose kinematical constraints to reduce uncertainties in both the experimental reconstruction 
of data and theoretical cross section of these QQqq states. This issue is the source of a 
recent discussion in the literature related with the secondary production of bb through gluon 
splitting 



5 Massless contribution to heavy quark sum rules 

Until present the evaluation of the perturbative two-point correlation function W ert (q 2 ) 
(in this Section we will denote the heavy quark currents correlator by U(q 2 )) has only been 
carried out completely, with massive quarks, up to 0(a 2 ) |J and the sum rules procedure, 
given by Eq. (H), has been termed consistent and effective in its task because the first branch 
point is set at the massive two-particle threshold. However Groote and Pivovarov have 
pointed out jlj that at 0(a 3 s ) there is a contribution to the correlator which contains a three- 
gluon massless intermediate state (see Fig. 0(a)). Its absorptive part starts at zero energy 
and, therefore, Eq. @ is no longer correct because on the r.h.s. there is a discontinuity 
starting at s = 0. Moreover those authors have also warned about the fact that, at this 
perturbative order, the massless intermediate state invalidates the definition of the moments 
Ai n for n > 4 because they become singular. Let us collect their reasoning. 

The perturbative contribution given by the diagram in Fig. Q(a) has been calculated at 
small q 2 (q 2 -C M 2 ) in Ref. [I]. In this limit the quark triangle loop can be integrated out and 
it ends up in the diagram in Fig. 0(b) generated by an induced effective current describing 
the interaction of the vector current with three gluons, 

J" = f (5 W + 14 W) , (12) 

with 

Oi U = d a bc G^ 13 Ga/3 , (13) 

where G% v is the gluon strength field tensor. The effective current in the QED case (G% u — ► 
F^jOig — > a e m,d a bc — » 1) can be easily identified from the lowest order Euler-Heisenberg 
Lagrangian (see Ref. |fj]). 

The correlator of the induced current (0) is then evaluated in the configuration space 
giving : 

<0|TJ„Cn) 4(0) |0> = - 2Q25 ^ 4M8 {^) 3 d abc d abc (d,d u - g^d 2 ) ± . (14) 
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Figure 2: (a) 0(a^) diagram contributing to the vacuum polarization function of the heavy 
quark current (the vertical dashed line indicates the massless cut), (b) "Effective" diagram 
obtained by integrating out the fermion loops. It also has the topological structure of the 
"reduced" diagram that determines the massless cut singularity. 



In momentum space we need to perform the Fourier transform of Eq. (0). Following the 
differential regularization procedure |J, which works directly in configuration space, the 
result for the vacuum polarization contribution of the diagram in Fig. ^](b) at small q 2 reads 



IW<7) 



17 



2916000vr 2 



dnhr.d. 



abcU'abc 



O 



3 / g 2 



V 



In 



+ o 



W 2 ^ 



(15) 



with n the renormalization point in this scheme, and d a b c d a bc = 40/3. 

As noticed by Groote and Pivovarov moments associated to the diagram in Fig. ^|(b) 
are not defined if n > 4. Indeed differentiating Eq. (|l^) four times, at q 2 ~ 0, we get: 




n(g 2 )| £ 



17 



218700tt 2 



4M 2 



In 



-q 2 



25 
12 



+ o 



M 10 



(16) 



whose real part clearly diverges if we set q 2 = 0. Larger n moments are also infrared divergent, 
and so the authors of Ref. [|IJ conclude that the standard sum rule analysis must limit the 
accuracy of theoretical calculations for the n > 4 moments to the 0(a 2 ) order of perturbation 
theory. This is, essentially, the conclusion of Ref. ||T| . 

An infrared safe redefinition of the moments, to cure the latter problem, has been provided 
in Ref. 0; it consists in evaluating the moments at an Euclidean point q 2 = —se, se > 0, 
thus avoiding the singular behaviour. This solution, as explained by the authors of that 
reference, is rather ill-conditioned from the phenomenological side though. Nevertheless the 
fault in Eq. (|8[) due to the massless threshold still represents a problem because even if, up 
to 0(al), we substitute the dispersion relation by 



W ert {q 2 



q 2 r ds hnngJ( S ) 



7T JAM 2 S 



+ 



)€ 



q 2 r°° ds Imn 39 (s) 
ir Jo s s — q 2 — ie 



(17) 
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(where Imn^(s) includes discontinuities starting at s — AM 2 ), the spectral function 
Imn 3s (s) associated to the cut in Fig. 0(a) would hardly be implemented phenomenolog- 
ically as gluons hadronize to both heavy and light quark pairs. We wish to provide a bypass 
to recover the balance between the right-hand and left-hand parts of Eq. flT7|). We will now 
see that if one does not insist in using full vacuum polarization for the sum rule analysis 
there is a way to overcome this infrared problem. 

In the heavy quark correlator the discontinuity across the three-gluon cut gives a contri- 
bution to the spectral function that is unequivocally real : 



i [U(s)] 3g = Imn 39 ( S ) = -j-J dR 3g {0\r\3g){3g\jl\0) 



from which the dispersive part can be evaluated independently of the QQ cuts. Accordingly 
we conclude that we can identify and isolate the troublesome massless cut contribution to 
the two-point function. Indeed Eqs. (|) and ([18]) justify our previous Eq. (|T/j). 

Let us go back then to Eq. ([17|) . All the difficulty with the phenomenological application 
of the sum rules is now the fact that the contribution from the three-gluon cut is contained 
in both sides of the equality. This intermediate state hadronizes completely into hadrons 
with a content of light and/or heavy quarks indistinctly. It is conspicuous that if we could 
disentangle the heavy quark hadronization, 3g — > QQ, we should include only this piece into 
the sum rule. Then the singularity at q 2 = would disappear because heavy quarks are 
produced starting at q 2 = AM 2 . However there is no way to sort out light and heavy quark 
production off three gluons and, therefore, if we extract this contribution from the heavy 
quark sum rules we are introducing an incertitude in the procedure because we make sure 
that there is no light quark hadronization but we miss the heavy quark production. It is easy 
to see that the induced error is small, due to the fact that three gluons hadronize mostly to 
light hadrons. On one side, in the very high energy region and following perturbative QCD 
with Np = 4, we have only a 1/4 = 25% probability of finding a specified pair of heavy 
quarks produced. And this is a generous upper limit because when we go down in energy, 
phase space restrictions severely reduce the counting of heavy quarks. Hence we estimate 
that excluding the three-gluon cut we introduce a tiny very few percent error in the sum 
rules procedure. 

Thus we propose an infrared safe definition of the moments by the trivial subtraction : 

fp-V) ^ firv, _ I r jf ton %M = tr jf Imn f (3) , ( i9) 

n Jo s s — q 2 — ie n Jam 2 s s — q 2 — ie 



u n ^ M.-ir, fa Im y . (20) 



Of course Eqs. (|19|) and (pip are meaningless unless we give a precise prescription about 
how to subtract the contribution of the massless cuts represented by ImH3 9 . Our previous 
discussion gives us the tool to proceed. Once the full 0(a>1) H pert (s) is calculated we can 
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extract the imaginary part starting at s = (which should go with a 9(s) function) for any 
value of s. It is clear that the 9(s) and 9(s — 4M 2 ) terms in the imaginary part of the vacuum 
polarization function correspond to three-gluon massless and to QQ cut graphs, respectively, 



and ImII 39 and Imll^ 4 are easy to distinguish, as Eq. (18) prevents the appearance of mixed 



9(s) ■ 9(s — AM 2 ) terms. Therefore we identify ImH3 g and we now plug it in the dispersion 
integral of the right-hand side of Eq. (|20|) and perform such integration. Divergences con- 
tained in both this integral and Ai n as q 2 — > will cancel with each other if the same infrared 
regularization is employed in the two quantities. The intuitive choice would be a low-energy 
cutoff s > 0, and Eq. fl2"0| ) would be more precisely written as: 



M n = lim 



1 I d Vrrperi/ 2m 1 f°° ds lmU 3g (s) 



U pert (q 



n\ \dq 2 I v ^ nq s ° ix Jo s (s + Sq) 



(21) 



where a vanishing term in the So — > + limit has been omitted. 

The evaluation of the M. n moments at q 2 = < AM 2 made sense because, up to 0(a 2 ), 
this point, being far away of the heavy quark production threshold, is unphysical and the 
moments are well defined through an analytic continuation from the high-energy region. 
However note that the absorptive three-gluon contribution starts at q 2 = and perturbative 
QCD becomes unreliable. This introduces a further new difficulty in evaluating Ai n moments 
at q 2 = 0, as we reach the physical non-perturbative region. Our definition of the moments, 
Ai n in Eq. (p0|) , skips this problem by fully eliminating the massless terms and, therefore, 
the final heavy quark sum rule will only involve physics at q 2 > AM 2 , apart from possible 
bound states. 

The general rule given above is valid for all orders of perturbation theory, but it strongly 
relies in our ability to extract the massless absorptive part from the full result of n(g 2 ) 
calculated at a definite order. Beyond 0(a 2 ) complete analytical results for the heavy quark 
correlator would be cumbersome and only numerical approaches may be at hand. In this 
sense, it would be convenient to have a method to calculate Imll^ only based on Feynman 
graphs. We have already sketched such a method in the discussion following Eq. (g) : we 
just need to sum up all the massless cut graphs to get Imn 3s , and then proceed with the 
dispersion integration that gives the associated dispersive part jlOH . For example, at 0(a^), 
the only massless absorptive part comes from the three-gluon cut in the diagram of Fig. 0(a); 
let us call M-sg the amplitude producing three gluons from the heavy quark current at lowest 
order (i.e. through the quark triangle loop in Fig. |3]). The massless contribution to the 
absorptive part of the correlator is then: 

Imn 3g (s) = ~ j dR 3g M% ■ Ml ga , (22) 



with the three-gluon phase space integral defined as 

1 1 7T 2 r 
3! (27r) 5 4Wo 



1 1 TV 2 f s fs-si 

dRsg = ^tttt^ y- / ds! / ds 2 , (23) 



o 
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Figure 3: Feynman diagram for the production of three gluons at 0{o? s ). 

in terms of the invariants s\ = (k\ + A^) 2 — (q — k^) 2 and S2 = (&2 + k^) 2 = (q — k\) 2 , and ki 
being the momenta of the gluons. The real part would be obtained by integrating Eq. ( |22|) : 

s r°° ds lmU 3g (s) _ -s r°° ds r' r>->* 

nJo ~s s + s -m^fJo WT7o)Jo dSl Jo ^ M ^' M ^ ^ 

which, in principle, could be performed also numerically. The nth-derivative of relation (p^) 
respect to sq, in the limit so — >• + , would give the infrared divergent contribution that should 
be subtracted from the full moments, as dictated by Eq. fl2"ip . 

Finally, we would like to mention that using the non-relativistic expansion of the heavy 
quark correlator in sum rules analyses does not avoid this infrared problem, at least formally. 
The 0{a Z g) diagram of Fig. [| will be highly suppressed in the velocity expansion, following 
the non-relativistic effective field theory approach, and therefore it is not relevant in the 
corresponding heavy quark currents correlator. However such two-point function cannot 
describe the QQ spectrum for energies far from threshold and even when higher n-moments, 
which strongly enhance the threshold, are used, perturbative QCD is needed in order to 
implement the remaining high-energy region; the diagram of Fig. ||] has to be accounted for 
to include properly this input, and its discontinuity at s = cannot be obviated. This point 
is more clearly seen by noticing that, besides the resummations in (a s /(3) performed in the 
non-relativistic correlator, one could improve such expansion by adding the terms needed to 
reproduce the exact 0{a z s ) result n(g 2 ). 

6 Conclusions 

Heavy quark sum rules, relying in global quark-hadron duality, are a compelling procedure 
to extract information on the theory from phenomenology. However, as higher perturbative 
order analyses are performed, the consistency of the method demands the inclusion of novel 
features. While at 0(a s ) the correlator of two heavy quark currents gives the full perturbative 
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information, at 0(a 2 s ) we have noticed that a heavy quark QQ pair radiated from light quarks 
in a correlator of light quark currents should be considered. At 0(a^) the complexities grow 
with the essential role of non-symmetric correlators. Closely related with this situation is the 
feature recently pointed out by Groote and Pivovarov on the uneasy problem arising from a 
massless three-gluon discontinuity in the heavy quark current correlator at 0(a^). 

We have shown that rigorous results of the general theory of singularities of perturba- 
tion theory provide all-important tools to analyse the new contributions. The inclusion or 
exclusion of specific discontinuities in the perturbative side is shown to be feasible and the 
decision involves a clear definition of the observable input on the phenomenological side of 
the sum rules. 

A solution for the problem pointed out by Groote and Pivovarov at 0(a^ ) has been given. 
We conclude that the appropriate procedure to obtain information about the heavy quark 
parameters should make use of the infrared safe corrected moments, defined in Eq. ([H]), that 
now indeed satisfy the modified sum rule : 

_ i ,oc Imnj^(s) 

M n = - ds (25) 

7i Jam 2 s n+1 

where the right-hand side can be extracted from the heavy quark production cross section 
<r(e + e~ — > [QQ]). The incertitude associated to heavy quark hadronization of the three-gluon 
should be taken into account but it is shown to be tiny. 

The analysis we have carried out is completely general, relying in the theory of singularities 
of perturbative theory amplitudes only, and provides a sharp tool for the future analysis of 
heavy quark sum rules. 
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